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Abstract. Recent and ongoing improvements to hydrodynamic treatments at RHIC are extending 
the physics reach of hydrodynamics, and improving the phenomenology. Here, the links between 
technological improvements and the extension of physics are emphasized. 

Hydrodynamic models represent the foundation for RHIC modeling. Using experimental data to address 
the equation of state, viscosity or dynamics of the novel matter created at RHIC would be impossible without 
hydrodynamic treatments. The last few years have seen a renaissance in the technical development of hydro- 
dynamical calculations. Improvements have included: solving fully three-dimensional equations, incorporating 
viscosity, coupling to dynamic mean fields. In this talk I will both describe each technical improvement and 
describe the physical implications and phenomenological consequences of the additional functionality. 



1 Solving fully three-dimensional equations 

Most of RHIC phenomenology has been confined to the central unit of rapidity. By assuming Bjorken-Hwa boost 
invariance, many hydrodynamic calculations have eliminated the need to model motion along the beam (z) axis 
by exploiting the independence of the equations for translations in the variable 77, defined by: 

z = TSmh(ri), i = rsinh(77). (1) 

After hydrodynamic equations are translated from x,y,z,t to x,y,ri,T, it is assumed that all quantities are 
independent of 77 which makes the equations effectively two-dimensional. In the boost-invariant ansatz the 
rapidity y of the matter is set equal to 77 which corresponds to acceleration-less matter, an assumption only 
warranted if there is a broad rapidity distribution. Otherwise, pressure gradients along the z direction would 
alter the evolution, and even affect the matter emitted at mid-rapidity. 

For observables outside central rapidity, the extension to three dimensions is not an improvement ~ it is a 
necessity, and I will not elaborate on the obvious benefits. Within the central unit of rapidity, the significance of 
performing three-dimensional calculations is subtle. First of all, the rapidity distribution at RHIC is consistent 
with thermal sources whose velocities are confined with ^ 1.5 units of rapidity of y = 0. Since a source moving 
with 1.5 units of rapidity can still emit a few percent of its particles with zero rapidity, the diligent modeler 
should modulate the source distribution accordingly. Crudely, this can be accomplished by introducing a weight 
Q-v /2A g^j-^y integral over source functions, which in the Hwa-Bjorken limit are uniformly distributed in rj. 

Modulations of the uniform Hwa-Bjorken source densities as described above neglect the effects of acceler- 
ation, as they continue to equate the coordinate 77 with the collective rapidity. Since longitudinal acceleration 
requires a longitudinal pressure gradient, acceleration vanishes if the 77 dependence is slow, or equivalently, if the 
reaction covers many units of rapidity. However, as the characteristic scale is between 1.5 and 2.0, acceleration 
cannot be easily dismissed. 

Several groups are now performing hydrodynamic calculations in three dimensions |l|2|3j . However, rather 
than using a full three-dimensional calculation to demonstrate the effects of longitudinal acceleration, I will 
employ simple one-dimensional solutions, where transverse degrees of freedom are neglected. In those solutions 
[4] one can compare the Hwa-Bjorken solution to one whose initial conditions, at some finite proper time tq, were 
consistent with having the collective rapidity equal 77. In the upper panel of Fig[Tl one can see that the width 
of the rapidity distribution broadens substantially due to the acceleration. This has important implications for 
models of stopping, which in almost all cases assume the final rapidity distribution has not evolved from the 
stopping stage. 
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Fig. 1. In the upper panel, the rapidity width of the source distribution is shown for a numerical solution to one- 
dimensional hydrodynamics (circles) which assumed that collective rapidity was equal to rj at the initial time of to = 0.5 
fm/c. It closely matches an analytic solution (solid line) that assumed complete stopping but matched the velocity 
gradient of the numerical solution at tq. This broadening is in stark contrast to the Bjorken-Hwa ansatz (dotted line) 
which assumes a fixed spread. For the same three cases, the inverse velocity gradient is shown in the lower panel. Whereas 
in the acceleration-less Hwa-Bjorken case the inverse velocity gradient equals t, it underestimates r by of the order 10% 
when acceleration is included. 

The inverse velocity gradient plays a critical role in two-particle correlation analyses, where the longitudinal 
size i?iong is identified with the inverse velocity gradient, 

-Rlong = {dv /dzy ^Vtherm- (2) 

Here, vthorm is the longitudinal thermal velocity, which can be extracted from blast-wave analyses of spectra [5]. 
Given that i?iong is directly taken from correlation analyses, one can extract the velocity gradient, which given 
the Hwa-Bjorken ansatz, t = [dv / dz)~^ , allows one to extract the breakup time. However, as can be seen in the 
lower panel of Fig. [TJ this equivalence between r and the inverse velocity gradient leads to an underestimate of 
the breakup lifetime on the order of 10%. Since modelers always had a difficult time explaining how the sources 
grew to their large observed transverse sizes within what appeared to be 10 fm/c or less, the extra 10% is rather 
welcome. 



2 Incorporating Viscosity 

Given the attention paid to viscosity in press releases announcing a "perfect fluid" at RHIC, it was surprising to 
see that none of the hydrodynamic treatments had incorporated viscous effects until the last year. Viscous shear 
increases the transverse pressure at the expense of the longitudinal pressure, thus giving transverse acceleration 
an early boost. This jump starts radial flow, and like longitudinal acceleration described above, helps explain 
correlation data. Although it is well understood that viscous shear can reduce elliptic flow, it is not clear that 
elliptic flow is reduced in all instances, and might even be increased by strong shear effects at very early times. 

It is worthwhile to first review the definitions and physics of viscosity. Rather than repeating the usual 
differential equations, I will write down the form of the stress energy tensor in the frame where the matter's 
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velocity is small, 



e (e + T^x)vx (e + Tyy)vy (e + T:,^)vz 

(e + Tyy)Vy 



T 











T 
^yy 














(3) 



where the coordinate axis have been chosen to diagonalize the stress energy tensor. It is always possible to 
diagonalize T"-^ by boosting and rotating the coordinate system. The above form is easily derived by applying 
a small additional boost to the diagonalized form and keeping only those terms linear in the velocity. In this 
coordinate system the equations of motion daT"i^ , 

(i9tV • v)e = -T^xdxVx - TyydyVy - T^z^^T^^, (4) 

dtVi = -^^diTii, (5) 

are physically transparent. Each element Ta plays the role of the pressure in the i^^ direction. Equation ^ 
expresses that the work done in expanding a small element in the i"^ direction is TadV, and in Eq. ([5]) the 
acceleration in the i direction is governed by the gradient of Ta . 

Viscosity concerns the physics of how Tij differs from the equilibrium value P{e)5ij. For Navier-Stokes, the 
difference is governed by the ansatz for the form of T*-' in the frame of the matter, 

Ty = P% - ry {d^vj + djv, - (2/3)V ■v)~BV-v, (6) 

which should be true for sufhciently small velocity gradients. Here, the coefficients r] and B are functions of 
the energy density and are referred to as the shear and bulk viscosities. Since the rate of change of the energy 
density is proportional to V • w, the bulk term describes how the effective pressure averaged over directions, 
T = (1/3) ^ Tii, departs from equilibrium due to a changing energy density. The strongest reason for T to depart 
from equilibrium is near Tc where the chiral field, and perhaps the gluon condensate, are rapidly changing with 
e [6l7j . Bulk viscosity will be discussed more heavily in the next section. The shear term describes the anisotropy 
in Tij due to the anisotropy of the velocity gradient, i.e., the traceless part of Tij. 

The Navier-Stokes equations suffer from a few shortcomings. First, numerical solutions include super-luminar 
modes [5]. Secondly, forcing the stress energy tensor to be determined solely by the velocity gradient, without 
any influence from prior history or initial conditions, limits the range of physics that can be addressed. For 
these reasons, a class of equations referred to as Israel-Stewart equations have been considered and are being 
implemented by several groups |9I10I11I8] . In these approaches the elements Ty are treated as dynamical objects 
that exponentially relax toward Navier-Stokes values. In addition to the two new parameters describing the 
relaxation times, r,j and r^, solutions can also differ by varying the choice of how one scales the relaxing 
variable. For instance, one can define the variables as: 



_ T,,-{l/i)5,,Y.J^^ 



(7) 



which then decay towards their Navier-Stokes values, 

(8) 

The usual choice for hjj and hs is that they are set to unity. Another choice might be to set h — P. 
The choices can yield quite different behavior for large departures from equilibrium, which can occur for rapid 
expansions. For instance, if a component Ta is half of the equilibrium value when a rapid expansion lowered the 
energy density to the point where the pressure was half of the original value in much less than the relaxation 
time, choosing h — 1 would lead to Ta falling to zero to maintain the offset with the equilibrium pressure. 
In contrast, if h were chosen to be equal to the pressure, the same rapid expansion would maintain the ratio 
Tii/P = 1/2. For a departure from equilibrium due to a slowly relaxing mean field, the first choice /is = 1 might 
be physically justified since freezing the field freezes the pressure contribution to the field. In contrast, when 
one considers the shear contribution from particles that cannot maintain equilibrium due to a low scattering 
rate, it would be the ratio of the anisotropy of Tij to the kinetic pressure that should freeze out. In fact, the 
choice oi hri = 1 could even result in negative components of Tij which should be unphysical for a dilute gas. 
This would suggest setting /i^ to the kinetic pressure, or perhaps to the enthalpy h. 
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In addition to the choice for h, one could also vary the target values for the offsets from the Navier-Stokes 
values assumed in Eq. ^ . Rather than allowing the target values to become arbitrarily strong for large velocity 
gradients, as is the case for the Navier-Stokes equations, one could choose a target that saturates for large 
velocity gradients, e.g., &'*^''8«'t = 5"^*- tanh(&*'*''s°*/6^'^*). For instance, if the system behaves like a free gas, one 
might enforce the pressure to remain positive, or if the matter is dominated by classical color electric fields, one 
might require the components of to be between — e and e. 

Since the components of are dynamical objects in Israel-Stewart calculations, they can also accommodate 
arbitrary initial conditions. For instance, if the energy density is initially dominated by large electric fields, 
the stress energy tensor approaches a form, Tzz — — e, Txx — Tyy — e. From this perspective, Israel-Stewart 
hydrodynamics is not so much a theory, but rather a surrogate model. Whereas for modest expansion rates Israel- 
Stewart theories reproduce Navier-Stokes hydrodynamics, for large expansion rates, the forms and parameters 
of the equations could be adjusted to mimic the behavior of a more sophisticated theory. For instance, if one 
is solving the equations of motion for the stress-energy tensor for Yang-Mills fields in a simple box geometry, 
one could tune the Israel-Stewart model to reproduce the evolution in the simple geometry. The Israel-Stewart 
model could then be tuned to more complicated geometries. Since the equation daT°'^ = is always true, 
and since varying the physical assumptions about the microscopic structure of the matter and equilibrium only 
affect the dynamics of Tij, any form of Israel- Stewart equations that correctly mimics the evolution of , will 
correctly replicate the evolution of the entire stress energy tensor and of the collective velocity. 

The importance of how one chooses to formulate the Israel-Stewart approach and its parameters becomes 
clear when one considers the strength of the Navier-Stokes corrections at early times. Even for the proposed 
minimal- viscosity limit, = s/47r, the Navier-Stokes correction to Tzz becomes 

s 4„ ^ / 4 \ 



T^-^ ^P-—ldzVz^P[l-—\, (9) 
47r 3 V STry ' ^ ' 

where the last step assumed an ideal gas where the entropy density is related to the pressure by P = 4sT, and r 
is the time. If the temperature is 300 MeV at r = 0.5 fm/c, then Tzz ~ P — — 1.75P, and the viscous correction 
is larger than the pressure itself. 

The principal manifestations of shear viscosity are at early times, where the large velocity gradients, dzVz = 
It, dxVx = dyVy ~ 0, is highly anisotropic. For free particles, the anisotropy at early times is bounded by the 
limit T^x = Tyy = e/2,Tzz = 0, whereas for longitudinal classical fields, the limit is T^x = Tyy = ^iTzz = ~e- 
For fields, the signs of the stress energy tensor can be understood by considering the expansion of the volume 
between two capacitor plates. Whereas pulling the plates apart requires work, expanding the area of the plates 
lowers the energy. Thus, during the first tq = 0.5 fm/c, the collective transverse velocities will expand by a 
characteristic amount 

where i? is a characteristic scale driving the transverse gradients. For off-center collisions R could be ^ 2 fm, v± 
could be between 0.8333 (particles) and 0.125 (longitudinal fields). It is clear that the initial transverse impulse 
in the pre-equilibrated stage is non-negligible and could provide 10-15% of the total collective flow generated 
during the expansion. The importance of this initial impulse in elliptic flow calculations remains to be analyzed 
in detail. 



3 Coupling to dynamic fields and bulk viscosity 

Viscous effects derive from a system's inability to maintain an equilibrium value, PSij, for the stress energy 
tensor at a fixed energy density, and a bulk viscosity ensues when '^.^Ta/S ^ P. For some cases, the bulk 
viscosity is manifestly zero. For instance, for a gas of free particles stress energy tensor becomes 

" E(^)v ' ^^^^ 

where i refers to the individual particles within the volume V. For a gas of ultrarelativistic particles Ta = e, 
regardless of the momentum distribution and Tu — e. The same can be said for the non-relativistic case once 
the particle number is fixed, and e refers to the kinetic energy density. Only a modest bulk viscosity will ensue if 
the distribution is semi-relativistic. Non-equilibrium chemistry can also lead to a bulk viscosity. For instance, in 
a cooling pion gas, the pion number would fall to maintain equilibrium as the temperature is near or below the 
pion mass |12|13j . Therefore, a rapid expansion, or equivalently a large V • v, will result in over-population of 
the pion number. As the extra pions reduce the net kinetic energy available at a fixed total energy, this can be 
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Fig. 2. The bulk viscosity due to non-equilibrium chiral fields are shown in the upper panels for two values of the 
parameter g from Eq. (|16p . For g = 3.4, the system is close to being a first-order transition and a changes rapidly as a 
function of energy density, and the bulk viscosity as calculated with Eq. (|15[) has a sharp peak. Assuming a simple one- 
dimensional Bjorken expansion provides the conditions for calculating the Navier-Stokes expectation for T = P — BV • v 
in the lower panel. The very sudden changes in T invalidate the Navier-Stokes approach and suggest that either an Israel- 
Stewart approach, or a simultaneous solution of hydrodynamic and field equations of motion, should be undertaken. 



considered as a viscous effect. However, if a model dynamically accounts for chemistry by assigning currents for 
non-equilibrated particle numbers, the non-equilibrated chemistry no longer contributes to the bulk viscosity. 

The most potentially important non-trivial sources of bulk viscosity derive from the vacuum condensates 
and fields which, near Tc, require a finite time to adjust to a rapidly changing energy density [5]. For a field 
obeying a differential equation, 

92^ + rat(/.-m2g((/,-^,q,ii)0 = O, (12) 
the drag term will pull the field away from the equilibrium value by an amount, 

r 

This result can be understood by considering the analogous problem of a harmonic oscillator moving through a 
medium that provides a drag force. Given that the equilibrium value of the field changes at a rate determined 
by the rate of change of the entropy density s — — sV • v, 

■ d(j)cquil _ /1 A\ 

(Pcquil = ' ^ ' 

one can identify the bulk viscosity as 



B = 



#oquil Sr dP 



ds mlg d4)eq 



(15) 



Since the effective mass will go to zero in the vicinity of a second-order phase transition, the viscosity will 
peak sharply at Tc [6] . Fig. [2] shows the bulk viscosity as a function of energy density for a quark gas in the 
linear sigma model, 

H = -Kv^a + ^ (^2 _ ^2 ^ ml/X'Y - h^a + i/quarks(m = ga), (16) 

After setting A — 1/40, the system has a first-order phase transition for g > 3.554, and a cross-over for lower g. 
The bulk viscosity, scaled by the pressure, is displayed in the upper panels of Fig. [5] for g = 3.4 and g = 3.2. As 
expected, the peak in the bulk viscosity is stronger for g near the critical value as that corresponds to a more 
rapid change in the a condensate as a function of energy density. In the lower panel, the effect on T is shown 
for the two values of g in a simple one-dimensional Hwa-Bjorken example. It is clear that viscous effects can be 
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large, and for the g — 3.4 case, so large that Navier-Stokes approaches are questionable. Solving Israel-Stewart 
equations should effectively smooth out the bumps in T, but as emphasized above, the form of the equations 
and the parameters should be carefully matched to reproduce the expected dynamics for the fields. 

Rather than solving Israel-Stewart equations, another option is to dynamically solve the equations of motion 
of the chiral field. This was done in [T3], though the equations ignored the drag term which from the consid- 
erations above should be important. The one crucial requirement for Israel-Stewart approaches to be valid is 
that equilibrium should be approached exponentially. In applying equations of motion for the field, the behavior 
is more of an oscillatory nature, although in the limit of strong damping the behavior becomes exponential. 
Thus, the appropriateness of Israel-Stewart approaches for modeling non-equilibrium behavior near Tc is not 
yet clear. However, explicity modeling the dynamics of the fields, alongside solving the hydrodynamic equations 
of motion, should always be tenable. 

We conclude with speculation regarding the effects of the non-equilibrium behavior near Tc. As can be seen 
from Fig. [21 the effect should be similar to what one would get from softening the equation of state near Tc, 
and perhaps even creating an unstable region with an effectively negative speed of sound. Lowering the speed 
of sound at the periphery of the collision volume has the effect of building up a density inversion, and thus a 
sudden dissolution of the system. This behavior should push the resulting Rout/Rside ratio from two-particle 
correlations downward, and perhaps in line with the low value observed experimentally. 
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